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SCALAR AND AXIAL MATRIX ELEMENTS OF THE NUCLEON: SEA QUARK CONTENT 
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Sea quark contributions to the scalar density and the axial current matrix elements of the nucleon are studied in 
lattice QCD with two flavours of dynamical Wilson fermions. The results are compared to trends in heavy quark 
mass expansions, and contrasted with the numbers obtained using dynamical staggered fermions. 



In the numerical study of many hadronic prop- 
erties, a comparison of the full QCD and the 
quenched approximation results shows that the 
differences are small (though the bare gauge cou- 
plings in the two cases are quite different); small 
enough that given the various sources of errors 
it is hard to differentiate between the two sets of 
numbers. Also the phenomcnological quark mod- 
els, which ignore the sea quarks altogether, do 
a good job of fitting many of the experimental 
results. All this suggests a welcome phenomeno- 
logical simplification of the theory (i.e. the dom- 
inant effect of the sea quarks is to generate con- 
stituent quarks as valence quarks with renor- 
malised properties), but one would like to un- 
derstand the dynamical reason behind it. 

To find unambiguous signatures of the sea 
quarks (or failures of the quenched approxima- 
tion), one has to look for instances where the sea 
quarks amount to more than mere renormalisa- 
tion of the gauge coupling. Good places to search 
are the properties where the quark model is not a 
good guide. Of particular interest are the quark 
bilinear matrix elements which naively vanish 
in the quark model, e.g. the strange quark ma- 
trix elements of the proton (p\sTs\p). We concen- 
trate on the scalar (T — 1) and the axial vector 
(r = 27^75) bilinears. Both are related to dy- 



namically broken classical symmetries of QCD, 
and hence have scope for a significant mixing 
with gluonic sector/sea quarks. Moreover, both 
correspond to Hermitian operators which are nu- 
merically easy to deal with, e.g. using the Hybrid 
Monte Carlo algorithm. (Complex fermion deter- 
minants arc difficult to handle as exemplified by 
the attempts to study QCD at finite chemical po- 
tential (F = 170) and the 9— dependence of the 
neutron electric dipole moment (F = 175).) 

The non-singlet quark bilinear matrix elements 
of the nucleon, conventionally expressed in terms 
of F and D couplings, are easily extracted from 
a direct evaluation of the 3— point correlators. 
On the contrary, correlators containing a purely 
gluonic intermediate state are quite noisy. Expe- 
rience has shown that instead of a direct evalu- 
ation of the 3— point correlators, it is easier to 
compute the sea quark matrix elements by mak- 
ing hadrons propagate through a background ex- 
ternal field (i.e. creating sea quarks with an extra 
source term Sr = J2x hrtp(x)Ti/}(x) added to the 
standard fermion action) and then evaluating nu- 
merical derivatives with respect to the external 
field strength Zip- While a singlet scalar density 
field is obtained just by choosing unequal valence 
and sea quark masses, configurations with a sin- 
glet axial current field have to be created afresh. 
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Pion-Nucleon Sigma Term 

The clear and interesting results concern the 
pion-nucleon sigma term (with m = m u = mj) : 

a vN = (N\m(uu + dd)\N) =m(dM N /dm) , (1) 

where the quark mass derivatives are to be eval- 
uated at fixed gauge coupling. The most recent 
analysis EJ of experimental data (though bet- 
ter data are definitely desirable) , gives u^n ~ 45 
MeV. Within the first order flavour SU(3) break- 
ing parametrisation, the valence quark compo- 
nent is only a$f = {3F S - D s ) w 26 MeV, leav- 
ing ample room for the sea quark component. 

We cannot yet directly calculate ct^n on the 
lattice, as it requires extrapolating the matrix el- 
ement to very small quark masses. Instead we get 
an indication of the importance of insertions on 
sea quark loops from the ratio of the full matrix 
element (valence plus sea) to its valence part. 

We generated 16 4 lattices with two flavours of 
dynamical Wilson fermions and periodic bound- 
ary conditions. These lattices were doubled in 
the time direction prior to calculation of hadron 
spectrum and matrix elements. We obtained re- 
sults § for /3 = 6g- 2 = 5.4, 5.5, 5.6, with the 
quark masses in the range m s < m q < 3m s . Our 
results for the ratio of the full matrix element 
to its valence part are shown in fig. 1. We see 
that the sea contribution is I — 2 times the va- 
lence part. This is in qualitative agreement with 
the experimental data. The overall magnitude of 
m&^i is systematically smaller than the experi- 
mental value, however, an effect possibly due to 
not having explored small enough quark masses. 

A similar analysis can be carried out for the 
staggered fermion results of the Columbia group 
||. Here the numerical derivatives for valence 
quark matrix elements involve large quark mass 
intervals, but the consistency between numbers 
extracted using (TV, SI) and (A, S) masses pro- 



Fig. 1. The results for sigma term matrix elements with 
two flavours of dynamical Wilson fermions. The curve is 
the leading term in lattice heavy quark mass expansion. 

vides a check on baryon masses being linear in 
quark masses over the interval. As table 1 shows, 
the effect of sea quarks is somewhat smaller 
in this case compared to Wilson fermions, but 
nonetheless significant. (Note that the effect is a 
factor of two smaller for heavy fermions.) 

The analogous scalar density matrix elements 
of the p and the A show similar factors between 
the valence and the full value, while the mag- 
nitude of the matrix elements is roughly pro- 
portional to the number of valence quarks. This 
suggests a model for constituent quarks in which 
the bare quarks are dressed strongly, and in a 
manner independent of the hadron that they are 
in. Indeed, it can be reasoned B that the large 
sea quark component we see is due to change in 
the overall scale of the theory; the sea quarks 
influence the /?— function through vacuum po- 
larisation. The really surprising feature of the 
lattice results then is (see fig.l) that even rel- 
atively heavy sea quarks (corresponding to the 
pseudoscalar meson mass up to 1 — 1.2 GeV) give 
a contribution comparable to the valence part. 
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Table 1 

The sigma term matrix elements with two flavours of dy- 
namical staggered fermions at /3 = 5.7. The baryon states 
used to calculate the matrix elements through numerical 
mass derivatives are indicated. 
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(A,S) 


0.010-0.015 




0.004-0.010 


9.0(3.) 


8.5(7) 


0.015-0.020 


19.2(2.8) 


0.015-0.070 


8.7(3) 


8.5(1.2) 


0.020-0.025 


11.6(2.3) 


0.020-0.070 


7.8(2) 


7.8(8) 






0.025-0.100 


7.4(1) 


7.5(3) 



Polarisation of Sea Quarks 

The EMC result on polarised muon-proton 
scattering is difficult to digest without a sizeable 
contribution from the sea quarks. Lattice mea- 
surement of sea quark and gluon components of 
the structure function g\ can be attempted in 
two ways. The simpler approach, possible even 
within the quenched approximation, is to replace 
the insertion on sea quark loops by an effective 
gluon operator Q. We have adopted the alter- 
native approach of directly determining the axial 
current coupling to sea quarks in the nucleon M . 

In the first trial run, we used two flavours 
of dynamical Wilson fermions on 8 4 lattices at 
(3 = 5.3 with k = 0.165,0.166. It is necessary 
to have the field strength Ka of the order of 
the quark mass to see any signal, and we used 
2nh,A — 0.005, 0.004. The pseudoscalar meson 
masses are about 0.9,0.8 GeV for these heavy 
quarks. We only find 2a bounds for the sea quark 
contribution to axial current matrix elements of 
the nucleon, \Sa\ < 0.05,0.08 per flavour. These 
are a factor of 2 to 4 below what is required to 
confirm the EMC result on the lattice. (The lat- 
tice one loop and continuum two loop renormal- 
isation constants for the singlet axial current are 
finite, and can be ignored at the present qualita- 
tive level.) This sharply contrasts with the scalar 



density case, where sea and valence quark com- 
ponents are comparable at similar quark masses. 
It is entirely possible that Sa moves away from 
zero at relatively small quark masses only. 

Lattice results for non-singlet scalar density 
and axial current matrix elements have turned 
out to be reasonable M . Heavy quark limits and 
constituent quark model expectations are easy 
to evaluate for these matrix elements. Putting 
everything together: (a) F$, mSs, Fa and Da 
are finite in the heavy quark limit and turn out 
to be smooth functions of the quark mass, while 
(b) mDs and Sa vanish in the heavy quark limit 
as mT 2 and depart significantly from zero only 
when the quark mass becomes of the order of 
h-QCD- This suggests that flavour SU(3) break- 
ing effects are probably small in the former case, 
though likely to be sizeable in the latter case. 
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